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Abstract: In a fractional Cauchy problem, the usual first order time derivative is replaced by 
a fractional derivative. This problem was first considered by Nigmatullin (1986), and Zaslavsky 
(1994) in M. d for modeling some physical phenomena. The fractional derivative models time 
delays in a diffusion process. We will give a survey of the recent results on the fractional Cauchy 
problem and its generalizations on bounded domains D C K d obtained in Meerschaert et al. 
(2009, 2010). We also study the solutions of fractional Cauchy problem where the first time 
derivative is replaced with an infinite sum of fractional derivatives. We point out connections 
to eigenvalue problems for the fractional time operators considered. The solutions to the 
eigenvalue problems are expressed by Mittag-Leffler functions and its generalized versions. The 
stochastic solution of the eigenvalue problems for the fractional derivatives are given by inverse 
subordinators. 
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1. INTRODUCTION 

A celebrated paper of Einstein (1906) established a mathe- 
matical link between random walks, the diffusion equation, 
and Brownian motion. The scaling limits of a simple ran- 
dom walk with mean zero, finite variance jumps yields a 
Brownian motion. The probability densities of the Brown- 
ian motion variables solve a diffusion equation, and hence 
we refer to the Brownian motion as the stochastic solution 
to the diffusion (heat) equation. The diffusion equation is 
the most familiar Cauchy problem. The general abstract 
Cauchy problem is dtu = Lu, where u{t) takes values in 
a Banach space and L is the generator of a continuous 
semigroup on that space, see Arendt et al. (2001). If L 
generates a Markov process, then we call this Markov pro- 
cess a stochastic solution to the Cauchy problem dtu — Lu, 
since its probability densities (or distributions) solve the 
Cauchy problem. This point of view has proven useful, for 
instance, in the modern theory of fractional calculus, since 
fractional derivatives are generators of certain (a-stable) 
stochastic processes, see Meerschaert and Shcffler (2004). 

Fractional derivatives are almost as old as their more 
familiar integer-order counterparts, see Miller and Ross 
(1993); Samko et al. (1993). Fractional diffusion equa- 
tions have recently been applied to problems in physics, 
finance, hydrology, and many other areas, see Gorenflo and 
Mainardi (2003); Kochubei (1990); Metzler and Klafter 
(2004); Scalas (2004). Fractional space derivatives are 
used to model anomalous diffusion or dispersion, where 
a particle plume spreads at a rate inconsistent with the 
classical Brownian motion model, and the plume may 
be asymmetric. When a fractional derivative replaces the 
second derivative in a diffusion or dispersion model, it 



leads to enhanced diffusion (also called superdiffusion) . 
Fractional time derivatives are connected with anoma- 
lous subdiffusion, where a cloud of particles spreads more 
slowly than a classical diffusion. Fractional Cauchy prob- 
lems replace the integer time derivative by its fractional 
counterpart: dfu = Lu. Here, d^g(t) indicates the Caputo 
fractional derivative in time, the inverse Laplace transform 
of s^g(s) — sP~ 1 g(0), where g(s) = / °° e~ st g(t)dt is the 
usual Laplace transform, see Caputo (1967). Nigmatullin 
(1986) gave a physical derivation of the fractional Cauchy 
problem, when L is the generator of some continuous 
Markov process {F(i)} started at x = 0. The mathe- 
matical study of fractional Cauchy problems was initiated 
by Kochubei (1989, 1990); Schneider and Wyss (1989). 
The existence and uniqueness of solutions was proved 
in Kochubei (1989, 1990). Fractional Cauchy problems 
were also invented independently by Zaslavsky (1994) as 
a model for Hamiltonian chaos. 

Stochastic solutions of fractional Cauchy problems are sub- 
ordinated processes. If X(t) is a stochastic solution to the 
Cauchy problem dtu = Au, then under certain technical 
conditions, the subordinate process X(E(t)) is a stochastic 
solution to the fractional Cauchy problem <9f u = Au, 
see Baeumer and Meerschaert (2001). Here, E[t) is the 
inverse or hitting time process to a stable subordinator 
D(t) with index (5 E (0,1). That is, E(t) = inf{x > 
: D(x) > t], and D{t) is a Levy process (continuous 
in probability with independent, stationary increments) 
whose smooth probability density fo(i){t) has Laplace 

transform e~ s>i = /d(i)(s), see Sato (1999). Just as Brow- 
nian motion is a scaling limit of a simple random walk, the 
stochastic solution to certain fractional Cauchy problems 



are scaling limits of continuous time random walks, in 
which the independent identically distributed (iid) jumps 
are separated by iid waiting times, sec Meerschaert et al. 
(2002). Fractional time derivatives arise from power law 
waiting times, where the probability of waiting longer than 
time t > falls off like for t large, see Meerschaert and 
Shcffler (2004). This is related to the fact that fractional 
derivatives are non-local operators defined by convolution 
with a power law, see Baeumer and Meerschaert (2001). 

In some applications, the waiting times between particle 
jumps evolve according to a more complicated process that 
cannot be adequately described by a single power law. 
Then, a waiting time model that is conditionally power law 
leads to a distributed-order fractional derivative in time, 
defined by integrating the fractional derivative of order 
(3 against the probability distribution of the power-law 
index, see Meerschaert and Shcffler (2006). The resulting 
distributed-order fractional Cauchy problem provides a 
more flexible model for anomalous sub-diffusion. The Levy 
measure of a stable subordinator with index (3 is integrated 
against the power law index distribution to define a subor- 
dinator W(t). Its inverse E(t) produces a stochastic solu- 
tion X{E{t)) of the distributed-order fractional Cauchy 
problem on M. d , when X(t) solves the original Cauchy 
problem du/dt — Lu. 

2. STOCHASTIC SOLUTION OF HEAT EQUATION 
ON BOUNDED DOMAINS 

Let I? be a bounded domain in M. d . We denote by 
C k (D),C k - a (D),C k (D) the space of k-times differentiable 
functions in D, the space of A;-times differential functions 
with fc-th derivative is Holder continuous of index a, and 
the space of functions that have all the derivatives up to 
order k extendable continuously up to the boundary dD 
of D, respectively. We refer to Meerschaert et al. (2009) 
for a detailed discussion of these spaces and concepts in 
this section. 

Suppose that every point of dD is regular for D c . The 
corresponding Markov process is a killed Brownian motion. 
We denote the eigenvalues and the eigenfunctions of the 
Laplacian A = Ylt=i d li by {/i„,^„}£L 1; where <j> n £ 
C°°(D). 

Remark 2.1. In particular, the initial function / regarded 
as an element of L 2 (D) can be represented as 



f{x) = y y\f{n)^ n {x)\ f(n)= / 4> n (x)f(x)dx. 

n=l Jd 



(1) 



The corresponding heat kernel is given by 

oo 

p D (t,x,y) = ^e _,tn Vn(aO0n(2/)- 



n=l 



The series converges absolutely and uniformly on [to, oo) x 
D x D for all to > 0. In this case, the semigroup given by 
T D (t)f(x) = E x [f(X(t))I(t < t d (X))] 

OO „ 

= y.e-'* nt <j>n(x)f(n)= / PD (t,x,y)f(y)dy 
n=i Jd 

(2) 

solves the heat equation in D with Dirichlet boundary 
conditions: 



d t u(t, x) = Au(t, x), xe D, t> 0, 
u(t, x) = 0, x £ dD, 
u(0,x) = f(x), xeD. 

3. FRACTIONAL CAUCHY PROBLEM 

Fractional derivatives in time are useful for physical mod- 
els that involve sticking or trapping, see Meerschaert et al. 
(2002). They are closely connected to random walk models 
with long waiting times between particle jumps, see Meer- 
schaert and Sheffler (2004). The fractional derivatives are 
essentially convolutions with a power law. Various forms of 
the fractional derivative can be defined, depending on the 
domain of the power law kernel, and the way boundary 
points are handled, see Miller and Ross (1993); Samko 
et al. (1993). The Caputo fractional derivative invented 
by Caputo (1967) is defined for < j3 < 1 as 



d^u(t, x) 



1 



r(i-/3) 7 

Its Laplace transform 



d r u{r, x) 



dr 



(3) 



f 

Jo 



e- st d?u(t, x) ds = s p u(s, x) - s p -^u{0, x) (4) 



incorporates the initial value in the same way as the first 
derivative. The Caputo derivative is useful for solving 
differential equations that involve a fractional time deriva- 
tive, see Gorenflo and Mainardi (2003); Podlubny (1999), 
because it naturally incorporates initial values. 

Let D £ R d be a bounded domain. In this section we will 
consider the fractional Cauchy problem: 

d?u(t,x) = Au(t,x), x£D,t>0; (5) 

u(t,x) = 0, x£dDt>0; 

u(0,x) = f(x), x £ D. 

To obtain a solution, let u(t, x) — G(t)F(x) be a solution 
of (5). Substituting in the PDE (5) leads to 

F(x)d?G(t) = G(t)AF(x) 
and now dividing both sides by G(t)F(x), we obtain 

afG(f) _ AFjx) 
G(t) F(x) A1 ' 

That is, 

df G{t) - —fiG(t), t > 0; (6) 
AF(x) = -(jlF(x) x£D; F{x) = 0, x £ dD. (7) 



Eigenvalue problem (7) is solved by an infinite sequence 
of pairs (/x ra , </>„), n > 1, where <p n is a sequence of 
functions that form a complete ortho normal set in L 2 (D), 
Mi < M2 < (J-2 < • ■ • , and fx n -> oo. 

Using the fi n determined by (7), we need to find a solution 
of (6) with ji = jj, n , which is the eigenvalue problem for 
the Caputo fractional derivative. 

We next consider the eigenvalue problem for the Caputo 
fractional derivative of order < j3 < 1 . 

Lemma 3.1. Let A > 0. The unique solution of the eigen- 
value problem 

d?G(t) = -\G{t), G(0) = 1. (8) 



is given by the Mittag-Lcfflcr function 

n=0 ^ 



pn) 



(9) 



For a detailed study of the Mittag-Lcffler type functions 
we refer the reader to the tutorial paper by Gorenflo and 
Mainardi (1997). 

Therefore the solution to (6) is given by 

G(t) = Go(n)M^(-/), 

where Go(n) = f(n) is selected to satisfy the initial 
condition /. Therefore using this lemma, we obtain a 
formal solution of the fractional Cauchy problem (5) as 



x ) = /( n ) M /3( _ M« 



(10) 



n=l 



Remark 3.2. The separation of variables technique works 
for a large class operators including uniformly elliptic 
operators L in divergence form, see Remark 4.4. For details 
see Meerschaert et al. (2009). Define a cube 

D = {x = (x\, X2, ■ ■ ■ , Xd) : < Xi < M for all 1 < i < d}. 

The functions 

d 

<f> n (x) = (2/M) d / 2 l[sm(nn i x i /M) 
»=i 

parametrised by the multi-index of positive integers n = 
{ni,ri2, ■ ■ ■ ,nd}, form a complete orthonormal set of 
eigenfunctions of the Laplacian, with Dirichlet boundary 
conditions with corresponding eigenvalues 

fi n = TT 2 M- 2 (n 2 1 + --- + n 2 d ). 

See, for example, Lemma 6.2.1 in Davies (1995). In this 
case the boundary of the cube domain is not smooth. 

3.1 Stochastic solution. 

Fractional time derivatives emerge in anomalous diffusion 
models, when particles wait a long time between jumps. 
In the standard model, called a continuous time random 
walk (CTRW), a particle waits a random time J n > 
and then takes a step of random size Y n . Suppose that 
the two sequences of i.i.d. random variables (J ra ) and 
(Y n ) are independent. The particle arrives at location 

X(n) = Yi H + Y n at time T(n) = Ji H + J n . 

Since N(t) = maxjji > : T(n) < t} is the number of 
jumps by time t > 0, the particle location at time t is 
X(N(t)). If EY n = and E[Y 2 ] < oo then, as the time 
scale c — >• oo, the random walk of particle jumps has a 
scaling limit c~ 1 / 2 X([ct\) => B(t), a standard Brownian 
motion. If P(J n > t) <~ ct~@ for some < (3 < 1 and c > 0, 
then the scaling limit c _1 / /3 T([ci]) =>• D{t) is a strictly 
increasing stable Levy process with index (3, sometimes 
called a stable subordinator. The jump times T(n) and 
the number of jumps N(t) are inverses {N(t) > n] = 
{T(n) < t}, and it follows that the scaling limits are also 
inverses, sec Meerschaert and Shcffler (2004, Theorem 3.2): 
c- f3 N{ct) E(t), where 

E(t) = inf{r : D(r) > t}, (11) 

so that {E(t) < t} = {D(t) > t}. A continuous mapping 
argument in Meerschaert and Shcffler (2004, Theorem 4.2) 
yields the CTRW scaling limit: Heuristically, since N(ct) s=y 



cPE{t), we have C -^ 2 X(N([ct})) w (c^)-^ 2 X(c^E(t)) w 
B{E{t)), a time-changed Brownian motion. The density 
u(t, x) of the process B(E(t)) solves a fractional Cauchy 
problem 

d^u(t, x) = d 2 u{t, x) 
where the order of the fractional derivative equals the 
index of the stable subordinator. Roughly speaking, if the 
probability of waiting longer than time t > between 
jumps falls off like t~P , then the limiting particle density 
solves a diffusion equation that involves a fractional time 
derivative of the same order (3. 

The Laplace transform of D(t) is given by 



-sD(t) 



e sx fD(t){x)dx = e 



The inverse to the stable subordinator E{t) of D(t) has 
density 

f E{ t){i) = §jPm) <i) = |(i - p(D(i) < t)) 



d rim 

diJo 



fD(i)(u)du 



(12) 



using the scaling property of the density /d(*)(0 = 
t~ 1/fj fD(i)(lt~ 1/f3 ), see Bcrtoin (1996). 

Using the representation (12) and taking Laplace trans- 
forms we can show that the unique solution of the eigen- 
value problem (8) is also given by 

/•OO 

G(t)= exp(-ZA)/ £(t) (Z)dZ = E(cxp(-A£(i))) (13) 
Jo 

see Meerschaert et al. (2002) for the details. 

Now using (2) and (13) we can express the solution to (5) 

as 

oo 

u(t,x) = f(n)Mp(-n n t )<f> n (x) 



: > f(n)[ / exp(-lfi n )f E(t) (l)dl}(t) n (x) 

oo 

[^2f(n)cxp(-ln n )(j) n (x)}f E{t) (l)dl \ > 



-f 

Jo 



n=l 



[T D (l)f(x)]f E(t) (l)dl 



= E x [f(B(E(t)))I(T M >E(t))] 

= E x [f(B(E(t)))I(T M {X(E)) > t)}. 

Remark 3.3. Meerschaert et al. (2009) established the 
conditions on the initial function / under which u(t,x) is 
a classical solution ( i.e., for each t > 0, u(t,x) e C 1 (D) n 
C 2 (D) and for each x £ D, u(t,x) £ C^O.oo)) of (5): 
that Af(x) has an eigenfunction expansion w.r.t. {4> n } 
that is absolutely and uniformly convergent. The analytic 
expression in (0, M) C K above is due to Agrawal (2002). 

4. DISTRIBUTED-ORDER FRACTIONAL CAUCHY 
PROBLEMS 

Let /jbea finite measure with supp^i C (0, 1). We consider 
the distributed order-time fractional derivative 



DM«(t,s):= / d?u(t,x)v(dp), v{dp) = T(l - /3)p,(d/3). 
Jo 

(15) 

To ensure that is well-defined, we impose the condi- 
tion 



< oo 



(16) 



as in Meerschaert and Shefflcr (2006, Eq. (3.3)). Since 
r(x) ~ 1/x, as x — > 0+, this ensures that v(df5) is a finite 
measure on (0, 1). 

4-1 Eigenvalue problem: solution with waiting time process 

Stochastic solution to the distributed-order fractional 
Cauchy problem is obtained by considering a a more flex- 
ible sequence of CTRW. At each scale c > 0, we are given 
i.i.d. waiting times (J n ) and i.i.d. jumps (Y£). Assume 
the waiting times and jumps form triangular arrays whose 
row sums converge in distribution. Letting X c {n) = Yf + 
■•• + Y° and T c (n) = Jf + ■ ■ ■ + J%, we require that 
X c (cu) A{t) and T c (cu) => W(t) as c ^ oo, where 
the limits A(t) and W(t) are independent Levy processes. 
Letting A t c = max{n > : T c (n) < t}, the CTRW 
scaling limit X C (N^) =>■ A(E") see, Meerschaert and Shcf- 
fler (2008, Theorem 2.1). A power-law mixture model for 
waiting times was proposed in Meerschaert and Sheffler 
(2006): Take an i.i.d. sequence of mixing variables 
with < Bi < 1 and assume P{J? > u\B t = 0} = c~ x u~& 
for u > c -1 /^, so that the waiting times are power laws 
conditional on the mixing variables. The waiting time 
process T c (cu) =>■ W(t) a nondecreasing Levy process, 
or subordinator, with E[e~ sW ^] = e -^ w ^ and Laplace 
exponent 



ipw{s) 
The Levy measure 

4>w{t, oo) = 



f 

Jo 



{e~ sx - l)(j) W (dx). 



(17) 



-Md/3), (18) 



/o jo r(i-/3) 

where fi is the distribution of the mixing variable, see 
Meerschaert and Sheffler (2006, Theorem 3.4 and Remark 
5.1). A computation in Meerschaert and Sheffler (2006, Eq. 
(3.18)) using f °°(l - e-^pt-P^dt = L(l - i3)s^ shows 
that 

lM*)= / ^L(l - finW) = f s^(d/3). (19) 
Jo Jo 
Then c~ 1 N£ => E v (t), the inverse subordinator, see Meer- 
schaert and Sheffler (2006, Theorem 3.10). The general 
infinitely divisible Levy process limit A(t) forms a strongly 
continuous convolution semigroup with generator L (e.g., 
see Arendt et al. (2001)) and the corresponding CTRW 
scaling limit A(E u (t)) is the stochastic solution to the 
distributed-order fractional Cauchy problem Meerschaert 
and Sheffler (2006, Eq. (5.12)) defined by 

{v) u{t,x) = Lu(t,x). (20) 
Since <j)w(0, oo) = oo in (18), Theorem 3.1 in Meerschaert 
and Sheffler (2008) implies that the inverse subordinator 
E v {t) = inf{x > : W{x) > t} (21) 

has density 



g(t,x) 



= [ 4>w{t-y,oo)P w{x) {dy). (22) 
Jo 



This same condition ensures also that E v (t) is almost 
surely continuous, since W(t) jumps in every interval, and 
hence is strictly increasing. Further, it follows from the 
definition (21) that E u (t) is monotone nondecreasing. 

We say that a function is a mild solution to a pseudo- 
differential equation if its transform solves the correspond- 
ing equation in transform space. The next Lemma follows 
easily by taking Laplace transforms. 

Lemma 4-1- (Meerschaert et al. (2010)). For any A > 0, 
h(t,X) = e- Xx g(t,x)dx = E[e- AB "(*)] is a mild solu- 
tion of 

B^h(t, A) = -Xh(t, A); h(0, A) = 1. (23) 

Kochubei (2009) considered the following: Let p(a) be a 
right continuous non-decreasing step function on (0,1). 
Assume that p has two sequences of jump points, /3 n 
and v n , n = 0, 1,2, where l3 n — > 0, v n — > 1, 
/?o = vq E (0,1). Suppose also that the sequence {/?«} 
is strictly decreasing and {v n } is strictly increasing. Let 
In = MM - pifin - 0)) and 7 2 = (p( Vn ) - p(v n - 0)), 
and define the distributed order differential operator 



u(t,x) = J2^^ 

n=0 



ln d^u(t,x) + J2lnd^u(t,x). (24) 

n=0 



Since p is a finite measure we have 

oo oo 
E 7n < °°> E In < °°- 



71=1 



n=l 



Here the corresponding subordinator is the sum of in- 
finitely many independent stable subordinators; 



wi = £(7i) 1/ ^r(i- / 8 n )w; 



n=0 



^( 7 2)V-r(i-^)VF- 



(25) 



71=1 



for independent stable subordinators Wf n , W" n for n = 

o,i,---. 

Lemma 4.2. Let EP(t) = inf{x > : W(x) > t}. 
Then h{t,p) — E(e _AtB "(^) is the classical solution to the 
eigenvalue problem 

B^h(t,p,) = -nh{t,n), h(0,p,) = 1. 

Using inverse Laplace transforms Kochubei (2009) estab- 
lished the following representation of h{t,p): 



h(t,p) 



a r r -i 

7f Jo 



where 
Hi(r) 

Mr) 



E 

n=0 



Mr) 

7 r l)r /3n sin(7r^„) + (7^K" sin(7ri/ n ) 



(26) 



E 



n=0 



y^r' 3 ™ cos(tt/3„) + ^r Vn cos(tw„) 



+ 



{OO 
n=0 



y^r' 3 ™ sin(7r/3„) + ^r Vn sin(7w„) 



(27) 



Let D C M d be a bounded domain with <9L> G C 1 '" for 
some < a < 1, and D x = (0, oo) x D. We will write 



u E C k (D) to meanjthat for each fixed t > 0,u(t,-) E 
C k (D), and u E C k (Ax,) to mean that u E C k (D 00 ) and 
is bounded. 

Define 

^a(-Doo) = {u : Ax, ->• R : Au(t, x) G C^)}; 
^i( J D 00 )='«A(£' 00 )n{u: |ft«(t,x)| < fc(t)ff(x), 
.9 G L°°(D), t > 0}, 
for some functions k and & satisfying the condition 
rl '* k(s)ds 



6(A) / / 
Jo Jo 

for t, A > and 



(*-*)/> 



< oo, 



k{t)Y^b{\ n )J{n)\^ n {x)\ < oo. 



(28) 



(29) 



n=l 



Theorem 1. Let / G C^-D) n C 2 (A for which the eigen- 
function expansion (of A/) with respect to the complete 
orthonormal basis {</>„ : n G N} converges uniformly and 
absolutely. Then the classical solution to the distributed- 
order fractional Cauchy problem 

B {p) u(t,x) = Au(t,x), xED,t>0; (30) 
u(t,x)=0, xEdD, t>0; 
it(0, x) = f(x), x E D, 
for u E % A (Ax>) n C&(Ax>) n C^-D), with the distributed 
order fractional derivative defined by (24), is given by 

u(t,x) =E x [f(B(Ef(t)))I(r D (B) > E"(t))] 
= E x [f(B(E?(t)))I(T D (B(EP))>t)} 

T D (l)f(x)g(t,l)dl 



(31) 



^2f(n)(t) n (x)h(t,ii n ). 



In this case, 6(A) = A, and k(t) is given by k(t) = Ci^ 0-1 , 

o < & < i. 

Proof. The proof is similar to the proof of Theorem 4.1 
in Meerschaert et al. (2010). We give the main parts of the 
proof here. 

Denote the Laplace transform t — > s of u(t, x) by 

/>oo 

u(s,x) = / e~ st u(t,x)dt. 
Jo 

Since we are working on a bounded domain, the Fourier 
transform methods in Meerschaert et al. (2002) arc not 
useful. Instead, we will employ Hilbert space methods. 
Hence, given a complete orthonormal basis {<p n (x)} on 
L 2 {D), we will call 



u(t, n) = / <fi n (x)u(t,x)dx; 
Jd 

u(s,n) = / 4> n (x) I e~ st u(t,x)dtdx 
Jd Jo 

= / cj) n (x)u(s,x)dx 
Jd 

/>oo 

= / e~ st u(t,x)dt (when Fubini Thm holds) 
Jo 

(32) 



respectively the </>„ and the ^„-Laplace transforms. Since 
{4>n\ is a complete orthonormal basis for L 2 (D), we can 
invert the </>„-transform to obtain 

M (*> x ) = ^2 u(t, n)i> n (x) 

n 

for any t > 0, where the above series converges in the L 2 
sense (e.g., see Royden (1968, Proposition 10.8.27)). 

Assume that u(t, x) solves (31). Using Green's second 
identity, we obtain 

4>n , du 



I [uA(j) n - <j) n Au}dx = 
Jd Jod 



86 



86 



ds = 0, 



since u\qd = = 4> n \dD, u E C l {D) by assumption, and 
<j> n G C 1 (D) by Gilbarg and Trudinger (2001, Theorem 
8.29). Hence, the <f> n -transform of Am is 



/ <fi n (x)Au(t, x)dx — — A„ / u(t, 
Jd Jd 



x)4> n (x)dx 



(33) 



= -X n u(t,n), 

as </>„ is the eigenfunction of the Laplacian corresponding 
to eigenvalue X n . 

The fact that the operator commutes with the <p n - 
transform follows from (28). 

Taking the </>„-transform of (31) we obtain that 

ID> (p) w(t,n) = -X n u(t,n). (34) 

From Lemma 4.2 we get the solution 

u(t,n) = f(n)h(t,X n ) = f(n)E(e' x " EP ^). 

Now inverting the ^-transform gives 

oo 

u{t,x) ='^2f(ri)(f> n (x)h(t,n n ). 
o 

The stochastic representation uses Lemma 4.2 and the 
stochastic representation of the killed semigroup of Brow- 
nian motion (2). 

We use the representation (26) to establish the fact that 
the solution is a classical solution. The details of the 
proof can be seen from the proof of the main results in 
Meerschaert et al. (2009, 2010). 

Remark 4-3. Suppose that n(d/3) — p(j3)dj3 1 the function 
p ^ r(l - P)p(p) is in C 1 [0, 1], supp(p) = [ft, fa] G (0, 1) 
and n{j3i) ^ 0. Then 



A 



_! _ tr $i(r) 
r e - ) , dr 



* 2 (r) 



h{t,X) =E[e- A£ " (t) 
where 

$i(r) = / r fi sin(/37r)r(l - P)p{fi)dp 
Jo 

$ 2 (r) = [ / cos(/37r)r(l - f3)p{fi)dp + A] 



(35) 



+ [f 1 rf j sm(/3n)r(l-l3)p(P)dp} 2 . 
Jo 

Suppose also that 



C(A),0i 



rPi 

p) = sin(/37r 

J Pa 



)r(l - P) P (f3)df3 > 0. (36) 



Then \d t h(t,X)\ < Xk(t), where 

k(t) = [CiPtoPup)*]- 1 ^! - + r(i - a,)**- 1 ]. 

(37) 

In this case, h(t,X) is a classical solution to (23). The 
representation (35) is due to Kochubei (2008), which 
follows by inverting the Laplace transform of (23). 

For u e H^iDoo) n Cbpoo) n C^-D) for k given by (37), 
Meerschaert et al. (2010) shows that the solution to (30), 
where 0( p ) replaced with the more general is a strong 
(classical) solution for / G C^-D) n C 2 (L») for which A/ 
has an absolutely and uniformly convergent eigenfunction 
expansion w.r.t {(j> n }- Nabcr (2004) studied distributed- 
order fractional Cauchy problem in D = (0, M) C M. 

Remark 4-4- The methods of this paper also apply to the 
Cauchy problems that are obtained by replacing Laplacian 
with uniformly elliptic operator in divergence form defined 
on C 2 functions by 



Lu 



E 



d(a ij (x)(du/dxi)) 
dxj 



(38) 



with a,ij(x) = dji(x) and, for some A > 0, 



"E^ 

i=l 



(39) 



*E y * ? - E a ij( x )yiVj ^ A ~ 

i=l ij — 1 

If A t is a solution to dX t = a(X t )dB t + b(X t )dt, X n = x , 
where <r is a d x <i matrix, and B t is a Brownian motion, 
then X t is associated with the operator L with a = aa T , 
see Chapters 1 and 5 in Bass (1998). Define the first exit 
time as td(X) — inf{t > : X t £ D}. The semigroup 
defined by T D (t)f(x) = E x [f(X t )I( TD (X)) > t)] has 
generator L with Dirichlet boundary conditions, which 
follows by an application of the Ito formula. 
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